Model of the quantized electromagnetic field in 
the presence of sources 

Jan Naudts and Wojciech De Roeck* 

Departement Natuurkunde, Universiteit Antwerpen UIA, 
Universiteitsplein 1, 2610 Antwerpen, Belgium; 
E-mail: Jan. Naudts® be and Wojciech.DeRoeck@ua.ac.be. 



Abstract 

We give a rigorous description of a model of the quantized electro- 
magnetic field interacting with quantized current fields. In the special 
case of classical currents our results agree with common knowledge 
about the problem. A toy model of a quantum current is studied as 
well. 



1 Introduction 

The traditional approach to quantum electrodynamics starts with the the- 
ory of free photons, respectively free electrons. Next, interactions are in- 
troduced. They are treated using scattering theory, making the assumption 
that particles are free in the large time limit. It is well known that the re- 
sulting theory suffers from intrinsic difficulties, like infrared and ultraviolet 
divergences. Here we show how to describe in a rigorous manner an electro- 
magnetic (e.m.) field interacting with quantized current fields. In principle, 
these fields might be generated by Dirac electron fields. However, we prefer 
to avoid the difficulties of fully quantized electrodynamics by delaying the 
rigorous formulation of electron fields to a forthcoming paper. The integra- 
tion of both models into a mathematically acceptable theory of QED remains 
an open problem. 
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Important for the rigorous formulation of the problem at hand is the 
choice of method for quantizing the e.m. field. As shown quite some time ago 
[1], the method of Fermi can be made rigorous. The essence of [1] has been 
picked up in [2] and was used to give an elegant formulation of quantized free 
e.m. fields using the formalism of covariance systems [3]. In particular, free- 
field operators are smeared out using test functions satisfying the continuity 
equation. By doing so they become gauge independent [4] and the difficult 
problem of gauge covariance of the quantized theory is avoided. 

Next a model is needed for the current fields. Given any such fields 
one can construct vector potential fields by integrating the currents with 
the Green's function of the d'Alembert equation. These vector potential 
fields can be described in a similar way as the free-field potentials. This is 
the basis for the present paper. The quantized e.m. fields and the currents 
are described simultaneously by a single covariance system. An essential 
property of such a covariance system is that only minimal information is 
needed about the properties of the two subsystems and of their interaction, 
because such information is encoded in the states of the system, which are 
determined by their correlation function. To illustrate this point we discuss 
states describing quantized e.m. fields interacting with classical currents and 
show that in this case we recover known results. But we discuss also states 
describing genuine interaction with arbitrary current fields. 

The structure of the paper is as follows. Sections 2, 3 and 4 serve as 
a technical introduction. Our basic ansatz is made in section 5. Section 6 
contains two propositions exploring the consequences of this ansatz. Sections 
7 and 8 show that our results are in agreement with results of the standard 
approach. In section 9 we introduce a simplified model of quantum currents. 
The physical interpretation of radiation fields of this model follows in section 
10. Finally, section 11 draws some conclusions. Two appendices contain 
some further technical matters. 

2 Test functions 

In photon theory a number of pitfalls have to be avoided and some choices 
have to be made. To begin with, it is well-known that the field operators 
Afj,(q) do only exist in a distributional sense. Therefore, smeared-out opera- 
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tors are defined by 

Mf)= I dqnq)A,(q). (1) 

The fn(q) are test functions reflecting the experimental inaccuracy to select 
a single point of Minkowski space. The Fourier transformed function 
satisfies 

f,(q) = (2vr)- 2 f dkf^e'^. (2) 

From the reality of f^(q) follows that f^{k) = f^(—k). 

Following [4, 2], we assume that the test functions satisfy the Fourier- 
transformed continuity equation 

WUW = 0- ( 3 ) 

A justification is given in Appendix A. This assumption is essential in the 
present formalism to control problems with gauge invariance of the theory. 

In what follows the space of test functions is denoted G. It consists 
of real functions f^(q) whose Fourier transform satisfies (3). For technical 
reasons, we assume that the Fourier transformed functions are continuous 
and vanish outside a bounded region. 



3 Classical wave functions 

It is tradition to introduce a complex Hilbert space of so-called classical wave 
functions of the photon. The space G is only a real space, because of the 
condition that fn(q) is real, while the classical wave functions form a complex 
pre-Hilbert space H. These are square integrable complex functions M (k), 
defined for k in M 3 , satisfying 

3 

|k|0 o (k)=^A; Q a (k). (4) 

a=l 

With each test function f^{q) corresponds a classical wave function M (k) by 
the relation 

M (k) = v / 27/ M (|k|,k), kGR 3 . (5) 
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The (degenerate) scalar product for classical wave functions is given by 
(m = -f dk-^-^(k)^(k). (6) 

JR3 Z\K\ 

Positivity of this scalar product follows because the classical wave functions 
satisfy (4). Indeed, one has 

(010) = - / dk-i-f|k| 2 |0 o (k)| 2 -|k| 2 ^|0 a (k)| : 
J* 3 2 l k l 3 V c^i 

= I dk ^lTl3 [^/3l k ! 2 - <Pa{k)4>p{k) 

a,p=l 

> 0. (7) 

The latter holds because the matrix with elements 5 a /?|k| 2 — k a kp is positive 
definite for all values of k. We conclude that the classical wave functions 
form a pre-Hilbert space H. 

4 Correlation function description of 
e.m. fields 

The smeared-out field operators, which will be constructed later on, satisfy 
the canonical commutation relations 

i(/),A(0)]_ = -2iIm<<^>, (8) 

where <fi and ip are the classical wave functions determined by / and g, 
respectively. The displacement operators are defined by 

W(f) = exp(ii(/)). (9) 

They satisfy the Weyl form of commutation relations 

W{f)W{g)=e nm ^W{f + g) (10) 

and generate an algebra which is not the algebra of canonical commutation 
relations [5] because the symplectic form 

/,<7-Im<^> (11) 
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is clearly degenerate. 

We need several distinct representations of this algebra. Free photon 
fields are described by field operators in Fock space. On the other hand, 
if the presence of an external current produces an infrared divergency, then 
a representation is needed which differs from the Fock representation. The 
obvious way to handle such a situation is by means of correlation functions. 
They determine the Hilbert space representation uniquely by means of the 
G.N.S. representation theorem. 

Instead of working with mathematical states of the C*-algebra of canon- 
ical commutation relations [5] we work with the correlation function formal- 
ism of [3]. The basic quantity is the two-point correlation function JF(/, (7) , 
defined for any pair of test functions / and g in G. In a Hilbert space 
representation with state vector Q it has the following meaning 

Hf,g) = {w(g)*n\w(fyn). (12) 

The scalar product between two vectors $ and ^ of the Hilbert space is 
denoted (\&| $) and is complex linear in $, anti-linear in ty. 

From (12) it is obvious that the correlation function J-{f,g) is nothing 
but the inner product between two coherent states, one with state vector 
W(f)*Q, the other with state vector W(g)*Q. The characterizing properties 
of correlation functions are in the present context 

• (normalization) JF(0,0) = 1; 

• (positivity) k \j\kJ r (fj 1 fk) > for all finite sequences of complex 
numbers Ai, • • • , A n and of elements /i, -f n of G. 

• (covariance) There exists a symplectic form a over G such that 

Hf + h,g + h) = e u ^-^(f,g) (13) 
holds for all /, g, h in G. 

The general definition of states of a covariance system contains also a 
requirement of continuity. However, the additive group of test functions 
G, + is equipped with the discrete topology. Hence continuity of the map 
/, g — > J-{fi g) is always satisfied. 

In particular, the correlation function for the vacuum state of the free 
photon field is given by [2] 

F(f,g) = e ^IV>} e -(i/2)<^} e -(i/2)<</>|V>) 
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= e ia(f,g) e -(l/2)s(f-gJ-g) 



(14) 



with 

<j(f,g) = lm(0|^) 

s(f,g) = Re (15) 

Here, and ip are the classical wave functions determined by / and g, re- 
spectively. Expression (14) satisfies all requirements for being a correlation 
function — see Appendix B. By the generalized GNS-theorem [3] there exists 
a projective representation W(f) of the group G in a Hilbert space H, and 
a vector f2 in H, such that (12) holds. 



5 Describing quantized currents 

The currents j^iq), given as operator- valued distributions in a Hilbert space 
H, can be used to construct vector potentials A? (q) by 

Al(q) = - J dq'A G (q'-q)Uq'), (16) 

where A G (g) is a Green's function of the d'Alembert equation, i.e. is a solu- 
tion of 

n q A G (q) = -5\q). (17) 

One can take Aa(q) equal to Feynman's propagator for massless bosons (47). 
The formal equation 

n q Al(q)=j,(q) (18) 

is satisfied by construction. After smearing out with test functions one ob- 
tains 

AKf) = -j d( l J dq'nq)A G (q' - q)j,(q f ). (19) 

Given these operators one can reconstruct the current fields by applying the 
d'Alembert operator. Indeed, introduce the notation 
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Then one has formally 



O a Ai(r a f) = J d q r(q-a)l(q) 

= 3(Taf). (21) 

Note that the free-field operators A(f) satisfy the homogeneous equation 

n a A(r a f) = 0. (22) 

As a consequence, two free-field operators A(f) and A(g) are equal if they 
determine the same classical wave function (modulo the null-space of H). 
This property is in general not true for the operators A*(f). Nevertheless 
one can use the formalism of covariance systems to describe these A*(f) in a 
similar way as for describing free field operators. Introduce the notation 

W*(/) = exp(iiJ(/)), (23) 

and assume that a correlation function J^(f,g) is given such that in the 
G.N.S. -represent at ion one has 

P(f,g) = {w i ( g yn\w*(fyn). (24) 

The current operators j(f) are fully specified by this correlation function. 
Examples of such functions follow below. 



6 Interacting fields 

Let us now construct an interacting field operator A 1 (q) which is the sum of 
the free field operator A^(q) and of the field A* (q) produced by the current. 
The latter two are described by the correlation functions J~(f,g), respec- 
tively JFJ (f,g). These have to be combined into a single correlation function 
F l {f,g), the G.N.S.-representation of which contains Weyl operators satis- 
fying 

W\f) = exp(ii I (/)) = exp(i(i(/) + AJ(/))). (25) 

An easy way to produce correlation functions with the desired properties 
starts from correlation functions JF X (/, /'; g, g') of the covariance system with 
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group G x G. By taking the diagonal of such a function one obtains a 
correlation function of the covariance system with group G 

F\f,g) = F x (fJ;g,g)- (26) 

The G.N.S.-representation induced by J~ l {f,g) can be obtained from that 
induced by F x (/, /; g, g). 

The simplest class of correlation functions of the product system consists 
of functions of the form 

F x (f,f;9,9') = e W (ia*(f,f;g,g')) 

x exp (-± 8 * (f -g,f- g'; f-g,f- </)) (27) 

with s x (f,f';g,g') a real inner product of G x G and with a x (f,f';g,g') a 
symplectic form of G x G such that 

(/, /'; g, 9') = * x (/, /'; g, 9') + i° x (/, /'; g, g') (28) 

defines a positive-definite sesquilinear form of G x G. These are analogues 
of the quasi- free states of [5]. Interaction between e.m. field and current is 
supposed to be such that 

F*(f,0;g,0) = F(f,g), 
.F x (0, /';(),</) = P(f',g')- (29) 

Let the G.N.S.-representation induced by J 7 * (/,/'; g, g') satisfy 

? x (/, /'; g, g') = (w*( g , <ft*n\ w*(f, /')*n) . (30) 

Introduce generators A l (f), A(f), and A j (/) by 

W*(fJ) = W I (f) = exp(iA I (f)) 
W*(f,0) = exp(ii(/)) 

^ x (0,/) = exp(^(/)) (31) 

By construction is ^4 ! (/) = A(f) +A*(f). The commutation relations for the 
operators W I (f) are 

W\f)W\g) = e^^W\f + g). (32) 
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The symplectic form appears in the r.h.s. of the commutation relations 

'A l (f),A\g)\_ = -2iaX(f,f;g,g) = -2ia\f,g) 
[A(f),A(g)]_ = -2^ x (/,0^,0) 

~A(f),A } (g)]_ = -2ta*(0,f;0,g) 

= -2m x (/,0;0,^). (33) 

7 Classical currents 

In the simplest case the currents j^(q) are multiples of the identity operator. 
Then the operators W^(f') and W ] (g') can be taken out of the inner product 
of (24). A suitable guess is therefore 

P(f,g) = exp(tA cl (g-f)) (34) 

with 

A c \f ) = - J dq Jdq' nq)A G (q' - q)j,(q'). (35) 

It is straightforward to verify that this function satisfies all requirements for 
being a correlation function. 

Eq. (34) allows very general classical potentials. Take e.g. the Coulomb 
potential 

At(q) = <WA (36) 

I M.I 

where c is a constant, and where q = (go 5 q)- Equations (35) are satisfied 
with external current 

j,(q) = ^,o47rc5 3 (q). (37) 

This means that a static charge of strength 47rc is located at the origin of 
space. The occurrence of a divergency in (36) does not produce any problem 
because it enters (34) in a form smeared out with test functions. 

Note that <J^(f,g) = 0. The commutation relations (33) suggest to take 

( T x (/,0;0,^) = ( T x (0 ) f;0,( ? ) = 0. (38) 
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It is still possible to include non-trivial correlations between the classical cur- 
rents and the free-field operators by means of the inner product s x (/, /'; g, g'). 
However, standard results about quantized e.m. fields in presence of classical 
currents are recovered when the choice 

(T x (f,f;g,g') = a(f,g) + A%-f) 

s x (f,f';g,g') = s(f, g ). (39) 

is made. 

All together, the correlation function of the classical current model reads 
F\f,9)=F(f,9)F(f,g) (40) 

with J-(f,g) given by (14) and fl(f,g) given by (34). Such a simple prod- 
uct form reflects the known fact that the classical current model does not 
contain any interactions. By this is meant that in a Heisenberg picture the 
Hamiltonian is the sum of two free parts, without additional interaction term. 

In the present model the interacting field operator equals the sum of 
the free-field operator and the classical potential generated by the external 
current 

A\f) = A(f) + A cl (f)L (41) 

This property is known in literature — see e.g. Eq. 2.63 of [6]. Still, many 
handbooks use the classical current model to illustrate the scattering ap- 
proach of QED, without mentioning (41). For sake of completeness we discuss 
some results of the scattering context in the next section. 

Note that the field operators A l (f) have some unusual properties. From 
(41) is clear that they satisfy the same canonical commutation relations as 
the free-field operators. But the representation depends intrinsically on the 
details of the external current. Indeed, a shift in spacetime may map a non- 
zero field operator onto zero. This means that the shifted representation is 
not unitary equivalent with the original representation. 

Let us analyze this point in somewhat more detail. A field operator A l (f ) 
vanishes if and only if A cX (f) = and (4>\ (f>) = 0, where is the classical wave 
function associated with /. Now, if the current j^(q) is not trivial, then there 
exists a test function / for which A cl (f) ^ and (0| 4>) = holds. If the 
current is localized in part of spacetime then shift the test function with a 
vector a so that r a f vanishes in that part of spacetime where the current does 
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not vanish. The result is that A cl (r a f) = 0. Because / and r a f determine 
the same classical wave function up to a phase factor one concludes that 
A\f) ± while A\r a f) = 0. 

8 Radiation fields 

Let us first verify what happens if A cl (f) is a solution of the homogeneous 
d'Alembert equation, i.e. the current vanishes. Then one can write 

A c \f) = ^dk^(^(k)^(k) + <k)0 M (k)) 

= -2Re(a|0), (42) 

with the classical wave function determined by the test functions / and 
with 

a M (k) = -l=^(|k|,k). (43) 

Hence the correlation function can be written as 

^(f,9) = ^(/^)exp(-2iRe(a|V-0)) 

= exp (ilm (0 + 2ia\ ip + 2ia)) exp(-(l/2)(0 -ip\(f)- ip)). 

(44) 

In this expression the Fourier transformed classical potential a, when multi- 
plied with 2i, behaves as a Fourier transformed test function. This relation 
is a duality between test functions and fields and has been studied in [2]. In 
particular, the definition of field operators A(f) can be extended to complex 
arguments. Hence one can write 

F\f,9) = F{f + nt- l A c \g + m~ l A d ) 

= (w{ g + in^A^yni W(f + in^A^yn^ 

= (Q(A cl )\W(g)W(fyQ(A cl )), (45) 

with 

n(A cl ) = w(i(27ry l A cl yn. (46) 
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To obtain the latter use that W(f + 2g) = W(g)W(f)W(g). This shows that 
in this case the correlation functions (40) are those of a coherent state, as 
expected from conventional photon theory. 

Next let us make a link with the scattering approach. Consider a classical 
field A c *(q) which vanishes for very negative times q « 0. This implies that 
jfi(q) = for very negative times. Assume that j^q) = holds also for very 
positive times go >> 0. Then A c ^(q) for g >> is the radiation field 
produced by currents Jn{q) which are only active during a finite interval of 
time go- This radiation field can be expressed in terms of the current using 
Feynman's propagator 

which is a Green's function of the d'Alembert equation. One finds 
K(l) = A T"(l) - J A l']M)D F (q - q') 

= A ^-<k*l ih ^r^~ m (48) 

where A l ^ om (q) is a solution of the homogeneous d'Alembert equation. Be- 
cause of the assumption that A^(q) = when g << 0, one must have 

= J^j^^f lkXqx Ukl go«0. (49) 

The latter equation can be written as 

4 om (g) = -J-Im / dk-^-e l l k l^e lk -%(|k|,k). (50) 
In J RS 2\k\ 

In this form the expression is valid for all g. Indeed, one checks immediately 
that this is a solution of the homogeneous d'Alembert equation. For go << 
one finds A^(q) = by construction, for g >> one has A c ^(q) = 2A^ om (q). 
Smearing out A^ om (q) with a test function / one obtains (see (42)) 

A ho m(/) = ^_ Im(0|a) (51) 

with the classical wave function determined by the test functions /, and 
with 

a M (k) = -^(]k|,k). (52) 



2tt 
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Hence, for test functions with support lying far in the future one has 

A c \f) = -lm( ( j ) \a). (53) 

71 

The correlation function (40) for a pair of such functions reads 

^(f,9) = ^(/^)exp((i/7r)Im^-^|a)). (54) 

Compare this with 

(W{g)*W{v)*tt\ W{f)*W{v)*Q) 
= exp(i Im (0 - V| 0) (W(v + g)*Sl\ W{v + f)*Q) 
= exp(ilm(0- ip\£))J 7 (v + f,v + g) 

= exp(2iIm<0-^|OW^), (55) 

where £ is the classical wave function corresponding with the test function 
v. The two expressions coincide provided there exists a test function v such 
that £ = — a/27r holds. If this is the case then the radiation field is described 
by the coherent state with wave vector W(v)*Q. This coincides with the 
standard result that, up to a phase factor, the S-matrix is a displacement 
operator, and the radiative field is a coherent state — see e.g. [7], section 
13. The problem is not the appearance of a complex test function, which has 
been explained above, but the possibility of an infrared divergency. Indeed, 
a test function v, such that £ = —ajlix holds, will not always exist. One 
can of course try to approximate a M (k) = j M (|k|, \i)/\/2ir by classical wave 
functions £ n . However, this will work only if 

(a\a) = _^jf^ dk^>(|k|,k)i(|k|,k) (56) 

is finite. This is precisely the condition for absence of infrared divergency 
found in [7], section 13, in case of an example. For recent progress on the 
infrared divergency problem in the context of Nelson's model see [8, 9] 



9 Quantum currents 



Let us now consider a simplified model of quantum currents. Start with 
creation and annihilation operators b* and b of a harmonic oscillator. They 
satisfy the canonical commutation relations 



b,b* 



(57) 
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Let be given complex functions a^(q) satisfying the continuity equation 

d^(q) = 0. (58) 
They are used to define currents j^(q) by the relation 

jM) = a M) h * + (59) 
The smeared-out potentials (19) become 

A>U) = yUy+W)b (60) 

with 

y(f) = " / d <? / d <?' n^a(q' ~ QhM')- (61) 
They satisfy commutation relations 

'A\f),Ai(g)]_ = -2ia\f ]g ) (62) 

with 

ai(f,g)= Im (fifty (gj) . (63) 

Let Q denote the ground state of the harmonic oscillator. It satisfies 
bfl = and determines the correlation function J^(f, g) via (24). One obtains 
by means of a standard calculation 

P(f,g) = e W (-iai(f,g))e W (-(1/2) j(f -g,f-g)) (64) 

with 

sKf, g) = \ (y(f)yjg) + v(g)vU)) • (65) 

In this state the quantum expectation of the current j^q) vanishes. The 
second moment equals 

(n\j^q)Uq')n) = ^Jq)a l/ (q f ). (66) 
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The interaction between photons and currents is modeled by assuming 
the existence of a real-linear function x(f) such that 



b,A(f) =x(f)i. 



(67) 



This implies 



v x U\0;0,g') 



i(/MV) 



= Im (x(f)y(g') 



The obvious choice of symplectic form a x (/, f';g, g') is then 
<r X (/, /';<?,</) = a(f,g)+^(f';g') 



+ lmx(f)y(g')+lmy(f')x(g)- 
The positivity requirement (28) suggest now to define 
s x (/, /';<?,<?') = s(f,g) + j(f,g') 



+ Rex(f)y(g')+Rey(f')x(g). 

Positivity is satisfied provided 

w)?<su\f) = m)- 

This implies the existence of functions and satisfying 

•(tf ) ,tf ) ) + «(/. Ca) , 

for which 



(68) 



(69) 



(70) 



(71) 



(72) 



(73) 



The symplectic form <? x (f, f; g, g') and the bilinear form s x (f,f';g,g') 
together determine the correlation function JF X (/, /'; g, g') via (27), and a 
corresponding state of the covariance system with group G x G. The diag- 
onal J~" l (f,g) describes a state of the quantized e.m. field interacting with a 
quantum current. It is the latter state which is analyzed below. 
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10 A quantum source of e.m. radiation 



First of all note that the state determined by J- l {f,g) from a classical point 
of view describes always a vacuum. Indeed, from 

^(/,0) = ^•( / ,0) e ( 1 /2)l-(/)l 2 e -(V2)|a ; (/)-y(/)| 2 (74) 
follows by expansion to first order in / that 

(A\f)) = 0. (75) 

The latter quantity is the classical part of the smeared-out e.m. vector po- 
tential. That it vanishes is in agreement with the pure quantum nature of 
the currents j^q) whose quantum expectation vanishes as well. 

Next consider the scattering situation with a current localized in space- 
time, i.e. cx^q) = outside some bounded region in the vicinity of the origin 
of spacetime. In addition, let the Green's function A G (g) in (61) be the 
retarded Green's function. Then y(f) = holds for all / with support in 
the far past. As a consequence, for /', g' with support in the far past is 
•T 7 * (/>/'> AS <?') — J~{f->g)- Hence, the state of the system in the past is the 
vacuum of the free e.m. field. 

The function y(f) does not vanish for all / with support in the future, 
even when a^(q) is again zero. In other words, the quantum current produces 
a radiation field. A first observation is that this radiation field cannot be 
coherent because the square of y(f) appears in (64). This is also obvious 
from (75). Nontrivial quantum fluctuations are present, as can be seen by 
expanding J- l {f,g) to first order in / and g 

(A\g)A\f)} = (m-^y(9)-yJf)<9) + yU)y(g)- (76) 

The first term in the r.h.s. describes the vacuum fluctuations. The last 
term describes field fluctuations which are identical to those of a classical 
radiation field. This leads to the remarkable observation that in this model 
quantum currents produce e.m. fields which propagate like classical radiation 
fields. They differ from them because their quantum expectation vanishes. 
Whether such quantum radiation fields exist in nature is not immediately 
clear. In more sophisticated models one can expect that these quantum 
radiation fields will survive in some sense. Indeed, any quantum current 
jfj,(q) can be decomposed into the sum of a classical current (j^(q)) and a 
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remainder with vanishing average. Now assume that the interacting field 
operators A l {q) satisfy the d'Alembert equation. Then, by linearity, the 
field operators decompose into the sum of a field produced by the classical 
current and a remainder without classical analogue, of the type found in the 
present model. 

A limitation of the model is that the quantum current is described by one 
single harmonic oscillator. As a consequence, the interacting field operators 
A 1 (f) satisfy unsatisfactory commutation relations (see (33) and (69)) 



A\f),A\g) = -2ia\f,g) (77) 



with 



Af,9) = a(f,g)-lmx(f)x(f) 



+ lm(x(f)+y(f))(x(g)+y(g)). (78) 

It is tempting to modify the model in such a way that the first two contri- 
butions to a l (f,g) cancel. One can hope that in such a modified model the 
energy density of the vacuum, which is infinite in absence of interactions, 
becomes finite in presence of quantum fields produced by quantum currents. 



11 Conclusions 

The present paper develops a rigorous theory for quantized e.m. fields in- 
teracting with given current fields. The formalism of covariance systems is 
used. The main tool in this approach is the vacuum to vacuum correlation 
function. It is not assumed a priori that the field operators are those of the 
non-interacting theory. Instead they are the sum of the free field operator 
and of a field operator which is the solution of the d'Alembert equation with 
the given quantum current as source term. 

The theory has been applied to two simple models, the first of which is 
the well-known model of quantized e.m. fields interacting with a classical cur- 
rent. We show that our approach agrees with standard results. In case the 
currents vanish outside a finite part of spacetime then our result describes 
radiation fields which correspond with coherent states, at least, if no infrared 
divergency occurs. The second model describes a current whose quantum ex- 
pectation vanishes. It produces a radiation field whose quantum expectation 
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vanishes as well. More sophisticated models are expected to produce similar 
results. 

Our main conclusion is that the algebra of field operators depends on 
details of the applied current. In this aspect our work goes beyond other 
approaches based on fixed algebraic structures and their representations. The 
emphasis on correlation functions solves the related technical problems in an 
elegant way. 

Field operators in the present paper are smeared out with test functions 
over spacetime, and not with test functions over 3-dimensional space, as 
is often done. We were not able to make a transition between these two 
approaches. In particular, we did not obtain a Heisenberg picture with a 
Hamiltonian, dependent on the external current, describing the time evolu- 
tion of field operators smeared out with test functions over 3-dimensional 
space. If it turns out that such description does not exist, then this is bad 
news for the standard approach, based on scattering theory, which takes the 
existence of an interaction picture for granted. 

The present work opens perspectives which may eventually lead to a 
rigorous formalism of quantum electrodynamics. The next step to take along 
the lines of the present paper is a rigorous description of a field of Dirac 
electrons interacting with a classical e.m. field. The algebra of Dirac currents 
is more complicated than what is supported by the present paper. The 
resulting technical problems have to be solved as well. 
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Appendix A 

Here we give a justification for assumption (3). 

Assume two vector potentials differ only by a gauge transformation 

A'Jq)=A ll (q) + d ll x(q), (Al) 
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with x(q) an arbitrary solution of the d'Alembert equation. Then the smear- 
ed-out fields satisfy 

A'(f)-A(f) = [ dqnq)d, X (q) 

= % [ &kx{-k)~r{k)k». (A2) 

Since the two vector potentials are physically equivalent, it should not be 
possible to distinguish them by means of the test function /. The condition 
A'(f) = A(f) implies then that (3) must hold for all wave vectors k satisfying 
= 0. 

Assume now that g^k) satisfies k^g a {k) whenever k v k u = 0, in such a 
way that the function 

X(k) = T^-k^(k). (A3) 

remains continuous. Then one can decompose g^k) into parts parallel and 
orthogonal to the wave vector 

Uk) = U(k) + k,x(k) (A4) 

where fn(k) satisfies (3) for all k. One verifies immediately that A(g) = 
A(f) holds for any vector potential A. Hence, we can always chose the test 
functions /, satisfying (3) for all k, as representative for a whole class of 
equivalent test functions g, satisfying (3) when k^k^ = 0, and such that (A3) 
remains continuous. 



Appendix B 

Here we show that the correlation function (14) satisfies the necessary condi- 
tions. The extension of the arguments to correlation function (40) is straight- 
forward. 

Normalization JF(0, 0) = 1 is clear. Positivity follows from 

EW(/ m ,/ n ) = ^A;A n e^l^e-^)^IW e -(V2)(/.l/.) 

mn mn 

= J2^n eiMfn) ( B1 ) 

mn 
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with 

/, m = A m e-^)<W/™>. ( B2 ) 

Note that the matrix with elements (f m \ f n ) is positive definite. Hence, 
positivity of (Bl) follows by means of Schur's lemma. 
Finally, covariance follows from 

jr(f + h ^ g + h ) = e i*n{f+h\g+h) e -{l/2){f-g\f-g) 

= e ilm{h\g) e ilm{f\h) :F ^^ g y (33) 
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